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ABSTRACT 

Context. The luminous infrared galaxy Zw 049.057 contains a compact obscured nucleus where a considerable amount of the galaxy’s 
luminosity is generated. This nucleus contains a dusty environment that is rich in molecular gas. One approach to probing this kind 
of environment and to revealing what is hidden behind the dust is to study the rotational lines of molecules that couple well with the 
infrared radiation emitted by the dust. 

Aims. We probe the physical conditions in the core of Zw 049.057 and establish the nature of its nuclear power source (starburst or 
active galactic nucleus). 

Methods. We observed Zw 049.057 with the Photodetector Array Camera and Spectrometer (PACS) and the Spectral and Photometric 
Imaging Receiver (SPIRE) onboard the Herschel Space Observatory in rotational lines of H 2 O, H^^O, OH, **OH, and [O I]. We 
modeled the unresolved core of the galaxy using a spherically symmetric radiative transfer code. To account for the different excitation 
requirements of the various molecular transitions, we use multiple components and different physical conditions. 

Results. We present the full high-resolution SPIRE FTS spectrum of Zw 049.057, along with relevant spectral scans in the PACS 
range. We find that a minimum of two different components (nuclear and extended) are required in order to account for the rich 
molecular line spectrum of Zw 049.057. The nuclear component has a radius of 10 - 30 pc, a very high infrared surface brightness 
(~- lO*'* Lq kpc““), warm dust (T^ > 100 K), and a very large H 2 column density (Nh, = 10^“* - 10^^ cm“^). The modeling also 
indicates high nuclear H 2 O (~ 5 x 10“^) and OH (~ 4 x 10“^) abundances relative to H 2 as well as a low *^0/*®0-ratio of 50 - 100. 
We also find a prominent infall signature in the [O I] line. We tentatively detect a 500 km s“* outflow in the H 2 O 3i3 ^2o2 line. 
Conclusions. The high surface brightness of the core indicates the presence of either a buried active galactic nucleus or a very dense 
nuclear starburst. The estimated column density towards the core of Zw 049.057 indicates that it is Compton-thick, making a buried 
X-ray source difficult to detect even in hard X-rays. We discuss the elevated H 2 O abundance in the nucleus in the context of warm 
grain and gas-phase chemistry. The H 2 O abundance is comparable to that of other compact (ultra-)luminous infrared galaxies such as 
NGC 4418 and Arp 220 - and also to hot cores in the Milky Way. The enhancement of '*0 is a possible indicator that the nucleus of 
Zw 049.057 is in a similar evolutionary stage as the nuclei of Arp 220 - and more advanced than NGC 4418. We discuss the origin 
of the extreme nuclear gas concentration and note that the infalling gas detected in [O I] implies that the gas reservoir in the central 
region of Zw 049.057 is being replenished. If confirmed, the H 2 O outflow suggests that the nucleus is in a stage of rapid evolution. 

Key words. ISM: molecules - Galaxies: ISM - Galaxies: individual: Zw 049.057 - Line: formation - Infrared: galaxies - Submil¬ 
limeter: galaxies 


1. Introduction 

Many luminous infrared galaxies (LIRGs) host compact ob¬ 
scured nuclei (CONs) where a bolometric luminosity Lboi > 10® 
Lq emerges from inside a core of diameter d < 100 pc that 
is obscured by dust corresponding to a high visual extinction 
Av > 1000 mag. The nature of the nuclear power source is thus 
hidden from examination with conventional methods like opti¬ 
cal and infrared (IR) lines and the Compton-thick shroud may 
render even X-rays unusable. It is important to determine if it 


* Herschel is an ESA space observatory with science instruments pro¬ 
vided by European-led Principal Investigator consortia and with impor¬ 
tant participation from NASA. 


is an accreting black hole or a compact starburst that powers 
the nuclear activity since this greatly affects our understanding 
of galaxy evolution. It has been suggested that a large portion of 
highly obscured Compton-thick active galactic nuclei (AGN) are 
missed by X-ray surveys, a problem that may be even worse for 
low-luminosity AGN (e.g.. lLusso et al.ll2013h . 

Furthermore, more tha n 50% of the star formation at high 
redsh ifts may be obscured (IChanman et akllZOO^ IWardlow et al.l 
I 2 OI ill . Obscured star formation can be linked to the assembly of 
stellar mass, with deeper potential wells in massive galaxies pro¬ 
viding dense, heavi ly obscured env ironments resulting in rapid 
star formation (e.g., iTharetaLlIZnnh . 
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The dust in obscured galaxies may have a high optical 
depth for a wide span of wavelengths. To study their nuclear 
structure, dynamics, and physical conditions, we thus need a 
tracer that can probe deep into the dust shroud. Molecules 
like water (H 2 O) and hydroxyl (OH) couple very well with 
the IR-field and can reach high abundances in warm em¬ 
bedded regions, making them ideal probes of physical con¬ 
ditions in dust-enshrouded galaxies as shown in previous 
studi es of Mrk 231, A r p 220 , and NGC 4418, among oth- 


20081 201QI 20121 l2014a 

Ibh. The OH molec 

:ule has also been 
1 ultraluminous in- 
Sturm et al.l201 ll 

proven to trace massive molecular outflows ii 
frared galaxies (ULIRGsl (iFischer et alJ2010 

Veilleux et al.ll2013t ISnoon et al. l2013l iGonzalez-Alfonso et al.l 

2014bh. Recent studies have found H 2 O in hi 

?h-redshift sources 

(Imoellizzeri et al. 2008; Omont et al.ll201 ll 

van der Werf et alJ 


the local Universe will thus aid any interpretation of such high- 
redshift observations. 

In the near Universe, a small sample of CONs have been 
identified with deep mi d-IR silicate absorp t ion and hot optically 
thick dust cores (e.g., [Aalto et^ l2012t ICostagl iola & Aaltol 


2 OIOI: C ostagliol a et alJ 2013 : Gonzalez-Alfonso et al. 12015 
Sakamoto et al.l ToiOi l2Q13ll . We have used the Photodetector 


Array Camera and Spectrometer (PACS; iPoditsch et al.ll201^ 
and the Spectral a nd Photometric Imaging Receiver (SPIRE; 
iGriffin et~^ 20 id) onboard the Herschel Space Observatory 
dPilbratt et al.l 1201*^ to observe H 2 O and OH in the CON Zw 
049.057, which is an OH megamaser dBaan et alJll987h and has 
mode rate IR luminosity (Lir 1.8 x lO" Lq H anderseT^ 
120031) . a very rich molecular spectrum in the f ar-IR and sub¬ 
millim eter, and a compact m olecular distribution (iPlanesas et al.l 
Il991h . Its NICMOS image (IScoville et al.ll20()0h reveals a nar¬ 
row dust featur e emerging along the minor axis from the ob¬ 
scured nucleus. IScoville et alJ (l2000l) also suggest that the fea¬ 
ture might be due to a nuclear absorbing cloud that blocks 
the light along the min or axis. Based on radio observations, 
iBaan & Klockneij (l2006l) classified Z w 049.057 as an AG N, but 
it is only weakly detected in X-rays (iLehmer et al.ll201^ . Mid- 
infrared ob servations with S'p/teer have revealed strong silicate 
absorption dPereira-Santaella et al.ll2010l) and far-IR ISO obser¬ 
vations indicated that the LIRGs Zw 049.057, NGC 4418, and 
IC 860 are all [C II] deficient dMalhotra et alJlT9^ . 

Our interest in this galaxy was originally sparked by the 
SPIRE FTS observations (see Sect. 12.41) . which revealed a rich 
molecular spectrum with H 2 O rotational lines that are compara¬ 
ble in strength to and, in some cases, stronger than the CO rota¬ 
tional lines. Modeling then raised the possibility that Zw 049.057 
has an extremely obscured nucleus with an embedded power 
source, possibly a young AGN. Very little is known about the 
obscured growth of AGNs, so that finding and identifying such 
objects is an important step in our understanding of galaxy evo¬ 
lution. 

The Herschel observations are described in Sect. ^ and our 
models are shown in Sect. [3] The model results are discussed in 
Sect.m and our main conclusions are summarized in Sect.|5] 


2. Observations and results 

Zw 049.057 was observed using the PACS and SPIRE spectrom¬ 
eters onboard the Herschel Space Observatory. The PACS obser¬ 


vation^ were conducted as part of the Hermolirg OT2 project 
(PI: Gonzalez-Alfonso) on 2012 July 20 in high spectral sam¬ 
pling range spectroscopy mode using first and second orders 
of the grating. The data reduction was done with the Herschel 
Interactive Processing Environment (HIPE; lOtd 1201 Oh version 
12.1.0 using the PACS background normalization pipeline for 
short range scans. At a distance of 56 Mpc, the far-IR nuclear 
emission of Zw 049.057 is spatially unresolved in the central 
9.1" (~ 2.5 kpc) spaxel of the 5x5 spaxel integral field ar¬ 
ray. Since the point spread function (PSF) is larger than the cen¬ 
tral spaxel, to obtain the full flux we extracted the central spec¬ 
trum using the point source correction task in HIPE 12.1.0. In 
each spectral range scan the continuum was then subtracted us¬ 
ing polynomials of third order or lower and the lines were fitted 
with Gaussian profiles. Line fluxes, measured and inferred in¬ 
trinsic Gaussian line widths, and continuum levels for the PACS 
observations are listed in Table[T] 

The corresponding SPIRE observation^ were conducted as 
part of the OT key program Hercules (PI: van der Werf) on 
2011 January 6 in high spectral resolution, single pointing, and 
sparse image sampling mode. A total of 108 repetitions (216 
FTS scans) were carried out, giving an on-source integration 
time of 14386 s. The data were processed using the standard 
point-source pipeline in HIPE 12.1.0. After subtracting the con¬ 
tinuum with a third order polynomial in each of the two bands, 
we measured the fluxes of the identified (> 3cr level) lines by 
simultaneously fitting the lines with line profiles consisting of a 
convolution of the FTS full resolution instrumental response (a 
sine function) with the best fit Gaussian line profile of the emis¬ 
sion from the galaxy. The full SPIRE FTS spectrum is presented 
in Sect. 12.41 Line fluxes and continuum levels for the H 2 O lines 
observed with SPIRE are listed in Table |2l In addition to the 
SPIRE observations, the Hercules program also included PACS 
observations of the [O I] 63 pm lin^l These observations were 
conducted on 2010 February 11. 

The spectral energy distribution (SED) of the source 
is shown in Fig. [1] including the SpitzerflRS spectrum 
( Armus e t al.l l2009l). and photometri c data points frorn IRAS 


( Sanders et alJ 200^ and SCUBA dPunne et al.l I 2 OO 0 I) . Note 


the strong silicate absorptio n at 9.7 and 18 pm reported by 
[Pereira-Santaella et al.l (l2010l) . 

Spectroscopic parameters for OH and H 2 O (both ortho and 
para) used for line identificati on and radiative tra nsfer model- 
ing were take n from the JPL (iPickett et al.lll998l) and CDMS 
(iMiiller et al.ll200ll l2005h catalogs. Energy level diagrams for 
the three species are shown in Fig. |2l where the transitions that 
were detected are indicated by blue and green arrows for PACS 
and SPIRE respectively. The detected transitions include H 2 O 
lines in absorption with lower level energies Eiower > 600 K 
(H 2 O 7 o 7 —> 6 i 6 and 7 i 7 —> 606 ), in emission with upper 
level energies also above 600 K (H 2 O 523 5 14 ),and OH 

lines with lower level energy Eiower > 400 K (OH ni /2 5 - 
i), indicating extreme excitat ion similar to that in Arp 220 
([Gonzalez-Alfonso et al.luM^ . Observations of highly excited 
H 2 O over the broad spectral range provided by the combina¬ 
tion of PACS and SPIRE give us an unique opportunity to 
constrain the parameters of the underlying continuum source 
with unprecedented accuracy. Note that the lines in the PACS 
range are all detected primarily in absorption while the lines 


' The observation identification numbers (OBSIDs) for the PACS ob¬ 
servations are 1342248365-1342248368 

2 The OB SID of the SPIRE observations is 1342203077 

3 The OBSID for the PACS observations of [O I] 63 pm is 1342190374 
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limeter wavelengths. Data from HerschelfPACS, /fericAeZ/SPIRE, and 
SpitzerflRS are shown. The IRAS data point s at 25, 60, and 100 ^ m as 
well as the upper limit at 12 pm come from lSan ders_et_^ 
(su b)millimeter points a t 850 and 1250 pm are from lPumie^ al.l ( I 2 OOOI 1 
and ICarico et alJ lll992h . respectively. The models discussed in Sect. [3 
are included where the blue and light blue curves represent the core and 
outer components, respectively. The green curve is an extra component 
to fit the continuum. Red is the total model. 


detected with SPIRE are mainly observed in emission, indi¬ 
cating that the rotational levels are pumped at far-IR wave- 
lengths and then relaxed by emissi on in the submillimeter 
(lGonzalez;Alfonso.e^^ 201(1 l2Q14ah . The analysis performed 
by iGonzalez-Alfonso et al.l (12014ah further indicates that the ab¬ 
sorptions are most efficiently produced close to the continuum 
source while the low-excitation (Euppei < 400 K) emission lines 
are formed in regions further away with less dust extinction, and 
as a result the combination of these wavelength ranges provides 
us with information on the source structure. 



excitation i n the lines indicates that they might be radia tively ex¬ 
cited (e.g.. iGonzalez-Alfonso et ^|2004 l2008l l2012ll and can 
be used to constrain the properties of the underlying continuum 
source. The spectral scans of the H 2 O lines are shown in Fig. [3] 
The absorption is not as strong as in NGC 4418, but is compa¬ 
rable to that of Arp 220 (iGonzalez-Alfonso et al.ll201^ . In ad¬ 
dition, nine transitions with upper level energies up to ~ 650 K 
were detected primarily in emission with SPIRE as summarized 
in Table |2] The spectral line energy distribution of the submil¬ 
limeter lines is shown in Fig. 0] As mentioned earlier, the emis¬ 
sion lines with upper level energies Eupper < 400 K have been 
shown to trace regions with low dust extinction, while the two 
detected emission lines with higher Eupper are found to form 
in very warm (> 80 K) regio ns with higher dust extinction 
(IGonzalez-Alfonso et al.l 1^14ah . Most of the lines are uncon¬ 
taminated by other species, but the H 2 O 3 12 —> 221 line is blended 
with the CO J - 10-9 transition and is not used in the analysis 
of the excitation. 


2.1.1. P-Cygni profile in H 2 O 3i3->2o2 

Although some lines appear shifted relative to the systemic ve¬ 
locity we find no systematic velocity shift, with some lines 
slightly shifted to the red and others to the blue. Also, within 2cr, 
all but two H 2 O lines peak at the galaxy redshift. The 3 i 3 —> 2 o 2 
line at 138 pm does however exhibit an apparent P-Cygni profile, 
possibly indicating an outflow. Signatures of this are not seen in 
any of the other PACS lines but out of our detected H 2 O absorp¬ 
tion lines, 3 i 3 ^ 2 o 2 is the one with the longest wavelength. It is 
thus the line that is least affected by dust extinction and it is pos¬ 
sible that the redshifted emission traces an outflow in the outer 
component on the far side of the galaxy that cannot be seen in the 
lines at shorter wavelengths, including the nearby 221 —»lio line 
at 108 pm, due to extinction. A possible explanation for the dif¬ 
ference in character of the two nearby lines, is that the 221 —> 110 
line, which is lower in energy, is also tracing the infall seen in 
the [O I] 63 pm line (see Sect. l2.31 l. in which case the redshifted 
absorption and emission could cancel. In fact, the profile of the 
221 —> lio line is slightly asymmetric, with less absorption on 
the redshifted side. The only H 2 O line detected in absorption in 
SPIRE (H 2 O 1 11 —> Ooo, see insert in Fig.[ 8 ]l tentatively exhibits a 
profile similar to the H 2 O 3 i 3 —) 2 o 2 line, with absorption at zero 
velocity and a redshifted emission component. 


2 . 1 . 2 . 

Two transitions in the isotopologue were also detected 

with SPIRE: 220 —> 2ii and 32 i ^ 3 i 2 with upper level energies 
of 194 and 303 K, respectively. While the 32i —> 3i2 line seems 
to be uncontaminated by other species, the 220 —> 2 ii is likely to 
be affected by the sine profile of the nearby H 2 O 422 —> 4 i 3 line- 


2.2. OH 


Fig. 2. Energy level diagrams of OH and H 2 O (ortho and para). Blue 
arrows indicate lines observed with PACS and green arrows denote lines 
observed with SPIRE. 

2.1. H 2 O 

Ten transitions of H 2 O with lower level energies up to ~ 650 K, 
were detected with PACS as summarized in Table [T] The high 


Five OH doublets with lower level energies up to > 400 K were 
detected with PACS, they are summarized in Table[T] The spec¬ 
tral scans of these lines are shown in Fig. |5] together with that 
of the undetected ni /2 3/2 1/2 doublet at 163 pm. Again, the 

absorptions are not as strong as in NGC 4418, but in the most 
excited lines they are almost comparable to those in Arp 220 
(IGonzalez-Alfonso et al.ll^l2ll . 

The doublet at 65 pm might have weak contamination by 
H 2 O 625 —> 5 i 4 in the blue A-component and the 84 pm dou- 
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Table 1. H 2 O, OH, and ‘'*OH lines detected with PACS 


Line 

■^rest 

(jjm) 

Ve"’" 
(km s“') 

(km s“') 

^^infeiTed 

(km S“^) 

Continuum*' 

(Jy) 

Peak line depth" 

(Jy) 

Flux" 

(Jy km s“') 

H 2 O 422 —>3i3 

57.636 

-25(17) 

319(43) 

247 

22.4 

1 . 8 ( 0 . 2 ) 

-613(108) 

H 2 O 524 —»4i3 

71.067 

35(39) 

275(98) 

221 

28.7 

0.7(0.2) 

-204(93) 

H 2 O 7i7 —>6o6 

71.540 

-12(65) 

292(147) 

243 

29.0 

0.5(0.2) 

-144(96) 

H 2 O 7o7 —»6i6 

71.947 

5(33) 

416(103) 

383 

29.3 

0.8(0.2) 

-358(112) 

H 2 O 321 —»2i2 

75.381 

20(4) 

235(9) 

178 

30.9 

6.4(0.2) 

-1595(76) 

H 2 O 423 —>3i2 

78.742 

-1(8) 

239(20) 

190 

31.2 

3.8(0.3) 

-975(105) 

H 2 O 6i6 —»5o5 

82.031 

-73(26) 

264(73) 

226 

31.0 

1.7(0.3) 

-464(158) 

H 2 O 6o6 —»5i5 

83.284 

-49(25) 

247(71) 

207 

31.2 

1.3(0.3) 

-339(123) 

H 2 O 221 —» lio 

108.073 

-19(10) 

366(26) 

194 

31.0 

4.0(0.2) 

-1541(145) 

H 2 O 313 —»2o2^ 

138.528 

-4(16) 

286(54) 

109 

27.8 

2.4(1.0) 

-716(325) 

H 2 O 313 —»2o2^ 

138.528 

266(173) 

508(219) 

434 

27.8 

-1.1 (0.4) 

606(334) 

OH n3/2 - n3/2 f ^ f 

65.132 

36(12) 

283(30) 

219 

26.6 

2.7(0.2) 

-817(111) 

OH n3/2 - n3/2 f - r 

65.279 

19(13) 

325(34) 

271 

26.6 

2.4(0.2) 

-833(114) 

OH ni/2 - ni/2 5 - f 

71.197« 

5(20) 

365(53) 

- 

28.7 

1.5(0.2) 

-596(111) 

OH ni/2 - n3/2 f - f 

79.118 

-22(7) 

148(18) 

- 

30.8 

3.1(0.3) 

-483(76) 

OH ni/2 - n3/2 - f 

79.181 

-18(7) 

125(19) 

- 

30.8 

2.9(0.4) 

-384(75) 

OH n3/2 - n3/2 f - f 

84.420 

12(4) 

227(9) 

185 

31.4 

7.2(0.2) 

-1742(89) 

OH n3/2 - n3/2 f - f ^ 

84.597 

27(4) 

221(9) 

178 

31.4 

7.3(0.2) 

-1712(90) 

OH n3/2 - n3/2 f ^ - f 

119.233 

37(5) 

253(11) 

- 

28.9 

6.8(0.3) 

-1842(102) 

OH n3/2 - n3/2 f ^ - r 

119.441 

13(4) 

236(10) 

- 

28.9 

6.7(0.2) 

-1695(93) 

'"OH n3/2 - n3/2 r - 

84.947 

8(9) 

99(23) 

- 

32.6 

1.8(0.4) 

-191(59) 

'"OH n3/2 - n3/2 r - 

85.123 

-39(20) 

201(48) 

- 

32.6 

1.3(0.3) 

-273(86) 

'«OHn3/2-n3/2i^-i^ 

120.171 

-36(33) 

254(82) 

- 

29.3 

0.8(0.2) 

-227(96) 


Notes. Values from Gaussian fits to the lines, numbers in parenthesis indicate Icr uncertainties from these fits. Velocity shift of line center 
relative to z = 0.012999. FWHM of lines. Inferred velocity width is based on the instrument resolution assuming a Gaussian profile; this is 
not listed for doublets that are not well separated. Value of the fitted baseline at the line center. ® The P-Cygni profile has been fitted using a 
combination of one absorption and one emission component. The two A-components are (nearly) blended into a single spectral feature. 


blet might be weakly contaminated by NH 3 ( 6 , 4)a (5,4)s and 2.4. Full SPIRE FTS spectrum 

( 6 , 5)a —> (5, 5)s in the red A-component. 


2.2.1. '®OH 

Two doublets of the isotopologue were also detected; 

ni/27/2 —> 5/2 at 85 pm and 113/25/2 —> 3/2 at 120 pm. 
The spectral scans of these lines are shown in Fig. |6l together 
with that of the undetected 113/2 9/2 —> 7/2 doublet at 65 pm. 
The blue A-component of the 120 pm doublet is strongly con¬ 
taminated by CH^ 3-2. Although the absorption in the main 
isotopologue generally seems to be weaker than observed in 
NGC 4418 by Gonzalez-Alfonso et al.l (l2()T3) . the '^OH absorp- 
tion is stronger in Zw 049.05 7. This is the same situation as 
iGonzalez-Alfonso et al.l (12012h found in Arp 220, likely indicat¬ 
ing enhancement of in Zw 049.057 as was inferred for Arp 
220 . 


2.3. [01] 63 pm 

The spectral scan of the [01] 63 pm line is shown in Fig.|7l note 
the prominent reversed P-Cygni profile, a typical signature of 
infalling gas. We do not attempt to model this line, but it provides 
important information about the kinematics on scales larger than 
the molecular components in our models. 


The full SPIRE FTS spectrum of Zw 049.057 is shown in Fig.O 
Because of the excellent match in the overlap region of the two 
spectrometer bands (~ 950 - 1020 GHz), they were simply av¬ 
eraged in this region. The spectrum shows at least 30 lines de¬ 
tected at a 3cr level or higher. Most of the CO ladde r is detected 
with nine lines from CO J ^ 4 - 3 to 7 = 12 - 11 (iGreve et al.l 
I 2 OI 4 I : l^senberg et al.ll2015h . In addition, nine rotational lines 
of H 2 O are detected together with two rotational lines of H^^O. 
The two [C i] fine structure lines and the [N ii] fine structure line 
as well as two rotational transitions of HCN are also detected. 
Finally, one rotational transition of NH 3 as well as one of HF, 
and at least three transitions of NH 2 are detected. There are hints 
of absorption accompanied with redshifted emission at the posi¬ 
tions of the two OH^ lines at 972 and 1033 GHz, but the noise 
in this region is too high for an unambiguous detection. Much 
stronger OH^ lines have been observed in emission in Mrk 231 
( van der Werf et al1l2010h and with P-Cygni profiles in Arp 220 
( Rangwala et al. 201 if) . The relative strength of the H 2 O lines 
with respect to the CO lines is greater in Zw 049.057 than in 
Mrk 231 and comparable to that seen in Arp 220. The lines are 
detected superposed on a continuum that drops towards the long 
wavelength side, and represents the Rayleigh-Jeans tail of the 
dust emission in Zw 049.057. 
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Table 2. HoO and Hi*0 lines detected with SPIRE. 


Line 

F 

^upper 

(K) 

drest 

(^m) 

Cont.“ 

(Jy) 

Flux" 

(Jy km s“') 


0.95 

H20 1n^0oo'^ 

53 

269.272 

6.9 

83.5(35.5) 


0.9 

H20 1ii^Ooo'^ 

53 

269.272 

6.9 

-146.8(35.5) 



H 2 O 2o2 ^ 1 11 

101 

303.456 

5.0 

904.1(53.3) 


1 

H 2 O 2ii ^2 o2 

137 

398.643 

2.2 

667.5(70.1) 

cd 


H 2 O 220 ^2ii 

196 

243.974 

9.0 

507.6(60.6) 

-i-j 

0 

0.95 

H 2 O 312 ^ 22 / 

249 

259.982 

7.7 

- 

Oh 

CO 


H 2 O 3i2^3o3 

249 

273.193 

6.7 

812.2(80.2) 

73 

0 ) 


H 2 O 321 ^3i2 

305 

257.795 

7.8 

1225.6(56.1) 



H 2 O 422 ^4i3 

454 

248.247 

8.6 

617.9(57.5) 

a 

0.9 

H 2 O 523 ^514 

642 

212.526 

12.6 

200.8(37.6) 

0 

0 

0.8 

H '«0 220 ^ 211 ^ 
H '«0 32i^3i2 

194 

303 

250.034 

263.738 

8.6 

7.4 

171.5(44.3) 

328.1(64.4) 

1 

a 

0 

0 

1.05 


Notes. Value of the fitted baseline at the line center. ® Values from 
fits of a convolution of a sine function and a Gaussian function to the 
lines, numbers in parenthesis indicate Icr uncertainties from these fits. 

The tentative P-Cygni profile in this line has been fitted using a com¬ 
bination of one absorption and one emission component. During the 
fitting, the frequencies for these components were held fixed at the po¬ 
sitions of maximum emission and absorption in the feature. Blended 
with CO y = 10-9, the total flux in the line is 1135.5(79.2) Jy km s“'. 

Likely to be affected by the sine profile of the nearby H 2 O 422 “^4i3 
line. 


1 

0,95 

1 

0.9 
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3 . Models 


As shown by the data presented in Sect. |2] the far-IR PACS 
spectra of Zw 049.057 are dominated by molecular absorp¬ 
tion, with only a few emission features, while the opposite is 
true for the submillimeter SPIRE spectrum. Because the ab¬ 
sorption and emission li nes of H 7 O require different phys- 
ical conditions to form (iGonzalez- Alfonso et alJ l2014ah . the 
galaxy must contain regions of different ISM parameters to 
produce these rich line spectra and their associated dust con- 
tinua. We attempt to model these regions using the smallest 
possible number of parameterized components. In doing so, we 
us e the spherically symmetric radiative transfer code described 
in I Gonzalez-Alfonso & Cernicharol (ll997ilT99^ . including line 
overlaps between between the A-components of the OH dou¬ 
blets in calculations for OH. Dust emission is simulated by us¬ 
ing a mixture of si licate and amo rpho us carbon grain s with opti- 
cal constants from lDrain^ d 19851) and iPreibisch et al.l (Il993h . the 
adopted mass absorption coefficient as a functio n of wavelength 
is described in iGonzalez- Alfonso et al.l (l2014ah. Rates for colli- 
sio nal excitation of H 2 O are taken from iDuternet et ^ (l2009h 
and Daniel et al.l (|201 ll) while the rates for OH were taken from 
lOffer et al.l (Il994l) . Note that the modeled regions are spatially 
unresolved by the Herschel beams. 

We find that at least two components are needed to fit the line 
absorption/emission and the continuum at wavelengths < 60 ^m. 
Depicted for the continuum in Fig. [T] a compact warm com¬ 
ponent (blue curve, hereafter Ccore) to account for the high- 
lying H 2 O lines and a larger cool component (light blue curve. 
Corner) to account for the low-lying OH lines as well as most 
of the flux in the submillimeter emission lines. The two com¬ 
ponents that account for all H2O and OH absorption/emission 
lines cannot, however, reproduce the observed continuum at 
d > 60 yum, so that an additional cold iji - 30 K) contin¬ 
uum component is added in Fig. [1] (green curve) to fit the whole 
SED. This extra component has no associated H2O/OH spec- 


Fig. 3. H 2 O lines observed in Zw 049.057 with PACS. The black his¬ 
tograms are the observed, continuum-normalized, spectra. Model re¬ 
sults from Sect. [ 3 ] are also included. The blue curve denotes the model 
for Ccore, which was constrained using the high-lying H 2 O lines detected 
with PACS. The light-blue curve denotes the model for Couter, which was 
constrained using the H 2 O lines detected with SPIRE. The red curve de¬ 
notes the sum of these two models. 



SPIRE, the black curve represents the data. For the lu —»Ooo line, only 
the absorption part is included. The line with upper energy level 249 K 
is 3 i2 —» 3o 3. Model predictions are included, the blue and light-blue 
curves show the contribution from Ccore and Comer, respectively. The red 
curve indicates the sum of the fluxes predicted by the two models. 


trum, which is expected given the relatively high required 
to excite these species. The parameters inferred for the two 
H 2 O/OH components are listed in Tables [3] and |4] In addition 
to these components, an even hotter and more compact compo¬ 
nent would be needed to fit the emission at wavelengths shorter 
than ~ 10 yum as was als o found for NGC 4418 and Arp 220 by 
IGonzalez-Alfonso et al.l (l2Q12l) . 
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Fig. 5. OH lines observed in Zw 049.057 with PACS. The black his¬ 
tograms are the observed, continuum-normalized, spectra. Model re¬ 
sults from Sect. [3 are also included. The blue curve denotes the model 
for Ccoro which was first constrained using the high-lying H 2 O lines 
detected with PACS, but with an added OH column that was varied in 
order to also fit the high-lying OH lines. The light-blue curve denotes 
the model for Comer, which was first constrained using the H 2 O lines 
detected with SPIRE, but with an added OH column that was varied in 
order to also account for the residual flux in the low-lying OH lines. The 
red curve denotes the sum of these two models. 


3. 1. Outline of the modeling 

In the following models we have assumed a standard gas-to- 
dust ratio of ~ 100 as foun d in the central regions of other 
LIRGs bv lWilson et al.l (l2008h . The high excitation in the molec¬ 
ular lines detected in absorption towards the nucleus cannot be 
accounted for by collisions alone, strongly suggesting that the 
dominant excitation mechanism is absorption of radiation emit¬ 
ted by warm dust a s in the nuclear regions of Mrk 231, Arp 
220 , a nd NGC 4418 (iGonzalez-Alfonso et al.ll2004l2008l l2()Ta 
I2OI2I) . Our models indeed indicate that collisional excitation in 
the core component described below is negligible in comparison 
with radiative excitation for reasonable densities and gas temper¬ 
atures (nHj ^ 5x10® cm“^, Tgas < 500 K), and thus our model re¬ 
sults are insensitive to Tgas. In the outer component, where the ra¬ 
diation field is weaker, collisional excitation is expected to have 
some impact on the low-lying submillimeter H 2 O lines (Sect. l3.3l 
below), but the excitation of the high-lying lines is dominated by 
radiative pumping. The excitation thus provides important clues 
about the properties of the continuum source and about the col¬ 
umn densities of H 2 O and OH. In all models the molecules are 
mixed with the dust, but as in Arp 220 (iGonzalez-Alfonso et al.l 
I 2 OI 2 I 1 the models indicate that most of the absorption is formed 
in the outer layers of the far-IR source where the continuum op¬ 
tical depth is < 1. The molecular column densities are therefore 
given for the outermost region in which T 50 = 1 , where T 50 is 
the dust opacity at 50 pm as measured from the observer. In case 
a model over predicts the flux density at 25 pm, it is attenuated 
by a foreground absorbing shell parameterized by its dust opac¬ 
ity at 25 pm. Evidence for such an abso rbing layer comes from 
the st rong silicate absorption at 9.7 pm (IPereira-Santaella et al.l 
I2OIOI) as well as from the absorption in the [O I] 63 pm line 
(see Fig. |7]). The total attenuated luminosity is however never 
allowed to ex ceed 10*^-^^ Lp, wh ich is the total IR luminosity 
of the galaxy (ISanders et alj|2003h . Details of the modeling are 
discussed in Sects. [3^ and l3.3l 



-1000 -500 0 500 1000 

Velocity (km/s) 


Fig. 6 . **OH lines observed in Zw 049.057 with PACS. The black his¬ 
tograms are the observed, continuum-normalized, spectra. Model re¬ 
sults from Sect. [3 are also included. The red curve denotes the model 
for Ccore. which was first constrained using the high-lying H 2 O lines de¬ 
tected with PACS, but with an added '*OH column that was varied in 
order to fit the **OH lines. 



Fig. 7. Line profile of the [O I] 63 pm line detected with PACS. The 
reversed P-Cygni profile in this line is a typical signature of infalling 
gas. 


In the models for the Ccore component, where collisional ex¬ 
citation is negligible, the line ratios depend on the following free 
model parameters; the dust temperature (Td), the column of dust 
(parameterized by the continuum optical depth at 100 pm, tioo), 
the column density of H 2 O per unit of T 50 (A^H 2 o/t 5 o), and the 
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N. Falstad et al.: Herschel spectroscopic observations of the compact obscured nucleus in Zw 049.057 
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Fig. 8. Full SPIRE FTS spectrum of Zw 049.057. Line identifications for > 3cr detections are given in red for CO lines, in blue for H 2 O, in magenta 
for H2*0, in black for HCN, and in green for the remaining lines. The insert shows a zoomed-in view of the H 2 O 1 11 —> Ooo line. 


velocity dispersion (AV, which is fixed, and derived from the 
observed line widths). In addition, absolute line and continuum 
fluxes depend on the size of the component (parameterized by 
the radius R, see Tables [3 and |4|l. As pointed out above, results 
of the models for the Couter component are also sensitive to the 
collisional excitation, and thus also depend on w h^ and Tpas- For 
these la tter parameters, we adopt the results by iMangum et al.l 
(I2013allil Tgas = 200 K and = 10^ cm-^). 

Our fitting strategy was the following: we developed a grid of 
models for H 2 O in both the Ccore and Couter components (some 
illustrative results are shown in Figs. l9] and [TOll by varying 
Tioo, and AIhjo/tso- Since both the Ccore and Couter components 
generally contribute to the flux of a given H 2 O line, we fitted all 
far-IR and submillimeter H 2 O lines simultaneously by combin¬ 
ing different parameters for Ccore and Couter- Once the best model 
fit for Ccore and Couter was obtained from all H 2 O lines, the same 
components were applied to OH by varying only the OH column 
in both components. 

3.2. The core component 

The best constraints on the dust temperature (Td) and column 
density (A^Hro) come from the high-lying lines of H 2 O, which 
are expected to be radiatively excited. In Fig. 0 we compare the 
observed and modeled fluxes of some of the highly excited lines, 
normalized to the flux in the 423 ^ 3 i 2 line. This line is chosen 


as the normalization because it is detected with high signal-to- 
noise and its lower level is high enough in energy for the line 
to remain relatively unaffected by extended low-excitation H 2 O. 
The modeled ratios are plotted as a function of Ta and A^Hjo/tso- 
We ran the models in Fig. |9]both with and without collisions in¬ 
cluded and found that the line ratios are unaffected by collisional 
excitation for nH 2 ^5x10^ cm~^ and Tgas < 500 K. For higher 
densities and gas temperatures the submillimeter H 2 O 1 n —»Ooo 
line will go into strong emission, which is not seen in the SPIRE 
spectrum. Line broadening is modeled by microturbulence with 
Vtur = 60 km s“5 

The relative fluxes of the high-lying lines can be well fitted 
with dust temperatures between 7^-90 and 130 K and column 
densities per T 50 between Ahjo/tso = 2x 10'* and 5 x lO'^ cm“^, 
although not all combinations of Ta and Nh^o/tso in this range 
produce a good fit. A lower dust temperature could also be used 
to fit the absorptions if a higher column density is used, but 
then the submillimeter emission lines 422 —> ‘^13 and 523 —> Su 
would be over predicted. In order to fit the absolute line fluxes 
the modeled region must have a radius between 32 and 12 pc. At 
these temperatures and radii, a foreground absorbing shell with 
an opacity between T 25 =2.1 and T 25 = 2.4 is required in order 
to fit the flux density at 25 pm. 

The dust opacity is best constrained by the high-lying sub¬ 
millimeter emission lines. With the column densities and tem¬ 
peratures required to fit the absorptions in the far-IR the emission 
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Fig. 9. Modeled H 2 O line ratios in Ccme as a function of dust temper¬ 
ature. Triangles, squares, pentagons, and circles show results for H 2 O 
columns per unit of T 50 of 2 x 10**, 1 x 10**, 5 x 10*^, and 2.5 x 10*’ 
cm“^, respectively. Dashed lines indicate the observed ratios and the 
dotted lines are their Icr uncertainties. 


line with the highest upper level energy, 523 —> 5 14 , would be over 
predicted for tiqo ^ 5 and under predicted for tiqo ^ 15. In or¬ 
der to ht the low-lying submillimeter lines of H 2 O an additional, 
less optically thick component is required, which is discussed in 
Sect. [33] 

To ht the high-lying OH lines we used the model parame¬ 
ters derived for H 2 O and varied the column density of OH to 
ht the observed line huxes. We hnd the best ht for OH/H 2 O ra¬ 
tios of 0.4 - 0.8, higher OH column densities over predict the 
ni /2 - ni /2 5 - I doublet at 71 pm while lower column densi¬ 
ties do not produce enough absorption in the n 3/2 - n 3/2 | - 5 
doublet at 65 pm. To ht the low-lying OH lines another compo¬ 
nent, which is discussed in Sect. 13.31 is required. The doublets 
of the isotopologue *^OH require high column densities in order 
to obtain a good ht. We hnd that an *®OH/'*OH ratio of 50 - 100 
is required to reproduce the absorption doublet at 120 pm while 
the one at 85 pm is slightly under predicted. 

The contribution of the core component to the total model 
is shown in blue in Figs. [T] and |3|6| In this reference model we 
have used - 120 K, A^Hjo/^so = 1 x 10'* cm“^, tiqo = 10, 
/? = 14 pc, 0 H/H 20 = 0.75, and '®OH/'*OH= 75. 

3.3. The outer component 

To adequately ht the emission in the submillimeter H 2 O lines 
and the absorption in the low-lying OH lines another, less ex¬ 
cited, component is required. This component also accounts 
for the remaining hux in the lower-lying H 2 O absorption lines. 


which are under predicted by the core component alone. In Fig. 
ITOlwe compare the observed and modeled huxes of some of the 
submillimeter lines, normalized to the Hux of the 2 o 2 ^ 1 11 line. 
The core component is expected to provide extra emission in 
the 2n —» 2 o 2 and 422 —> 4 i 3 lines. Some of the submillimeter 
lines, for example H 2 O 422 ^ 4 i 3 , are however still not repro¬ 
duced, probably indicating a transition region between the two 
components where the emission in these lines is still significant. 
In this component we hnd that the ratios of the lowest lying lines 
are significantly altered by collisional excitation, in the models 
shown here we have adopted a gas te mperature - 2 00 K and 
a density hhj = 10^ cm * as derived bv iMangum et al.l(l2013allbh 
based on observations of formaldehyde and ammonia. A poten¬ 
tial caveat is that they did not consider IR pumping when esti¬ 
mating their temperatures. The main effect of increasing the den¬ 
sity is to strengthen the low-lying H 2 O In ^ Ooo and 2 o 2 —»In 
lines at the expense of the H 2 O 32 i —> 3 i 2 line. With a higher den¬ 
sity the two lower lines would be over predicted while a lower 
density produces too much emission in the H 2 O 32 i —> 3 i 2 line. 

In order to reproduce the line ratios, we hnd that a dust tem¬ 
perature between Ta = 45 and 60 K and column densities per T 50 
between A^H 2 o/t 50 = 2 x 10 *^ and 8 x 10 *^ cm“^ are required. 
The line ratios are also sensitive to the dust opacity, we hnd that 
the ratios can be reproduced with tiqo from 0.5 to 2. In order to 
reproduce the absolute Huxes the component must have a radius 
of 90 to 50 pc. The line widths are still best reproduced using 
Vtur = 60 km s'*. 

The submillimeter H^^O 32 i —» 3 i 2 line is mainly formed in 
the outer component, requiring an abundance ratio *®0/**0< 50 
in water. With this ratio for OH, the '*OH 120 pm doublet would 
be overestimated by a factor ~ 2, but the **OH 85 pm doublet 
would be better reproduced than in our reference model (Fig. | 6 ]l. 
We thus favor a range *®0/**0= 50 - 100, with enhanced '*0 
abundance in both the core and outer components. 

The OH lines are again htted with the model derived for the 
H 2 O excitation. In this component a OH/H 2 O ratio of 0.1 to 0.4 
is required to ht the low-lying OH lines. With higher OH column 
densities the undetected doublet at 163 pm is over predicted. 

The contribution of the outer component to the total model 
is shown in light blue in Figs. [T] and |3|6| In this reference model 
we have used - 55 K, A^Hjo/'J'so = 4 x 10*^ cm“^, tiqo = 1, 
R = 60 pc, and OH/H 20 = 0 . 2 . 

3.4. Model uncertainties 

We have ignored all inner structure in the modeled components 
by adopting uniform densities and dust temperatures throughout 
the models. The approximation of spherical symmetry is also 
clearly inaccurate, the NICMOS image of Zw 049.057, for ex- 
ampl e, exhibits a linear str ucture extending for hundreds of par¬ 
secs (IScoville et al.l l2()00h . For the absorption lines, geometry 
should however be secondary as the important property is the 
projected surface of the underlying optically thick far-IR source. 
Our models should nevertheless be considered as very simplihed 
versions of the real structures. 

The dust temperature in Ccore is well constrained as high tem¬ 
peratures are required for the radiatively excited high-lying lines 
to form. In Comer we expect that collisional excitation is also im¬ 
portant and Td is here less certain than in Ccore- Due to dust ex¬ 
tinction, the absorption lines are sensitive to the molecular col¬ 
umn densities only in the outer layers of the core where Tgo < 1 . 
In these layers they are however fairly certain and present robust 
lower limits to the total column densities of the source. The con- 


Article number, page 8 of 1 141 


















































N. Falstad et al.: Herschel spectroscopic observations of the compact obscured nucleus in Zw 049.057 


Table 3. Parameters of the continuum models. 


C“ 

Radius* 

Td* 

Tioo* 

A’h/ 


U 

T25,fg/ 

Tatt« 


(pc) 

(K) 

at 100 yum 

(cm"’) 

(10^ Mo) 

(Lo) 


(Lo) 

c 

^core 

32-12 

90-130 

5-15 

(3 - 10) x 10’“* 

0.3-7 

(0.7- 1.2) x 10“ 

2.1 -2.4 

(2-6)x 10*** 

^^outer 

90-50 

45-60 

0.5-2 

(0.3- 1.3)x 10’“* 

0.5-7 

(4-9)x 10'° 

0 

(4-9)x 10*** 


Notes. Component. ® Independent parameter. Column density of H 2 , calculated assuming a mass-absorption coefficient of 44.4 cm^ g“* 
at 100 yum and a gas-to-dust mass ratio of 100. Estimated mass, assuming spherical symmetry. Unattenuated luminosity of the component. 
® Foreground opacity at 25 yum. Attenuated luminosity of the component. 

Table 4. Derived H 2 O column densities and abundances, and column density ratios 


C** 

A’h2o/T50*''* 

(cm"’) 

,, d 

XH 2 O 

OH/H 2 O*'’ 

16Q^18Q/j,e 

c 

'-'Core 

(0.5 - 2) X 10*" 

(2.5 - 10) X 10"* 

0.4 - 0.8 

50 - 100 

^outer 

(2-8)x 10*’ 

(1 -4)x 10"* 

0.1 -0.4 

50 - 100 


Notes. Component. ® Independent parameter. Column densities per unit of dust opacity at 50 /rm, T 50 . Estimated H 2 O abundance relative 
to H nuclei. Column density ratios. 



30 40 50 60 70 

Td (K) 

Fig. 10. Modeled H 2 O line ratios in Coutcr as a function of dust temper¬ 
ature. Triangles, squares, pentagons, and circles show results for H 2 O 
columns per unit T 50 of 8 x 10 *’, 4x 10 *’, 2 x 10 *’, and 1 x 10 *’, respec¬ 
tively. In all models Tjoo = 1. Dashed lines indicate the observed ratios 
and the dotted lines are their Icr uncertainties. The contribution from 
Ccore is not included in this figure, hut it is expected to slightly raise the 
2ii —> 2 o 2 / 2 o 2 —> 111 and 422 4 i 3 / 2 o 2 —> In ratios. The overall best fit, 
including the core contribution, is shown in Eig.|4] 


is propagated to the molecular abundances, which of course de¬ 
pend on the column of H 2 . Abundance ratios between different 
species are more certain as the various abundances are calculated 
in the same way. 


4. Discussion 


4.1. The dust and gas components ofZw 049.057 

As described in Sect. |3] above we can fit the spectral line data 
with two components whose properties are summarized in Ta¬ 
bles [ 3 ] and 0] For comparison, Table |5] contains a summary of 
physical conditions in some other galaxies as well as in some 
Galactic regions. 

Compared to giant molecular clouds (GMCs) in the Galaxy, 
Ccore is comparable in size but significantly hotter as well as 
more massive and luminous (IScoville & Goodl|l989t) . Although 
cooler and less dense, it shows strong s imilarity to the compact 
hot co re found in the LIRG NGC 4418 (iGonzalez-Alfonso et al.l 
l2012h . The physical conditions are also very similar to those in 
the western nu cleus of the ULIRG Arp 220, w hich is a few times 
as big as Ccore llGonzalez-Alfonso et al.ll2()T^ . 

The outer component is larger than a typical GMC and is 
(again) significantly more massive and luminous, its dust tem¬ 
perature is however co mparable to the peak te mperature found 
in H Il-region GMCs dScoville & Goodll989h . An interesting 
comparison is with the Galactic Center - where a mass of a 
few times 10’ Mp of molecular gas is d i stributed on a scale of 
200 - 500 pc ( Scoville & Sanders 198^ ICox & Laureiis 1989: 


Morris (fe SerabytiT 1996t iDahmen et al.l Il998t iFerri^re et all 


2007h . Thus the nuclear region of Zw 049.057 appears to be dra¬ 


matically different from the center of our own galaxy, with an 
order of magnitude more molecular gas concentrated on scales 
comparable to a large GMC. 


4.2. Abundances and column densities toward Zw 049.057 


version of tiqo to a column density of H 2 depends on the adopted 
gas-to-dust ratio as well as the dust properties. This uncertainty 


All abundances derived for Zw 049.057 are summarized in Table 
0] For comparison, a summary of abundance estimates in other 
galaxies and selected Galactic regions is tabulated in Table |5] 
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Table 5. Physical conditions in other objects for comparison 


Object 

Class 

Nh, 

(cm-2) 

Amo 

Aoh 

R 

(pc) 

Ld 

(K) 

M 

(10^ Mo) 

Ref. 

APM 08279+5255 

HyLIRG 

- 

6 X 10-'“ 

- 

100 

- 

- 

1,2 

Mrk 231 warm 

ULIRG 

- 

8.0 X 10-'^ 

- 

120 

95 

5.9 

3 

Mrk 23 1 extended 

ULIRG 

- 

- 

- 

610 

41 

77 

3 

Arp 220west 

ULIRG 

(3 - 6) X 10^5 

(0.4 - 3) X 10-5 

(0.16 - 3) X 10-5 

75-40 

90-130 

55-30 

4 

Arp 220east 

ULIRG 

1.3 X 10^5 

1.6 X 10-5 

~ 1.6 X 10-5 

54 

87 

13 

4 

Arp 220extended 

ULIRG 

8.0 X 10^3 

8 X 10-* 

~ 2.4 X 10-5 

325 

90-40 

30 

4 

NGC 441 Score 

LIRG 

(0.5- 1.0) X 10^5 

(1 -3) X 10-5 

(0.3 - 1.8) X 10-5 

10 

140-150 

0.16-0.33 

4 

NGC4418warm 

LIRG 

2.4 X 10^4 

(1 -5) X 10-5 

(0.4 - 2) X 10-5 

15.5 

110 

0.2 

4 

NGC 441 Sextended 

LIRG 

(1.0- 1.3) X 10^3 

4 X 10-* 

(2-4)x 10-5 

100 

90-30 

0.18 

4 

Zw 049.057core'' 

LIRG 

(3 - 10) X 10^4 

1 

0 

X 

0 

1 

a\ 

(1 -8) X 10-5 

32- 12 

90-130 

0.3-7 

5 

Zw 049.057outer'' 

LIRG 

(0.3- 1.3) X 10^4 

(1 -4)x 10-5 

(0.1 - 1.6) X 10-5 

90-50 

45-60 

0.5-7 

5 

Object 

Class 

Nh, 

(cm“^) 

AhjO 

Aoh 

R 

(AU) 

Ta 

(K) 

M 

(10« Mo) 

Ref. 

Sgr B2 

Denser regions 

- 

few X 10-' 

- 

- 

- 

- 

6 

Sgr B2 

Envelope 

- 

(1 -2)x 10-5 

(2 - 5) X 10-5 

- 

- 

- 

6,7 

Orion KL 

Hot core 

- 

1.0 X 10-5 

- 

- 

- 

- 

8 

Orion KL 

Plateau 

- 

7.4 X 10-5 

(0.5 - 1.0) X 10-5 

- 

- 

- 

8,9 

Orion-IRc2 

Shocked low-velocity gas 

- 

(2 - 5) X 10-4 

(6.2- 15.6) X 10-5 

- 

- 

- 

10 

NGC 6334 I 

Hot core 

- 

~ 10-5 

- 

- 

- 

- 

11 

NGC 1333 IRAS2A 

Hot core 

- 

> 2 X 10-5 

- 

~ 100 

- 

- 

12 

W43-MM1 

Hot core 

- 

1.4 X 10-4 

- 

~ 3500 

- 

- 

13 

L1157-B1 

Warm shocked region 

1.2 X 10^“ 

(0.7 - 2) X 10-5 

- 

1250 

- 

- 

14 

L1157-B1 

Hot shocked region 

3.3 X lO^o 

(1.2 - 3.6) X 10-4 

- 

250 - 625 

- 

- 

14 


Notes. The value from ivan der Werf et alJ ll201 ih has been used. 


References, til [ B radford et al.l 1 I 20 IT ): 121 [ v an der Werf et al.l ( 201lh: t3~) iGonzalez-Alfonso et al.l (l2010h: (41 iGonz alez-Alfonso et al.l (1201 ^; (5) This v york; (6) ICernic haro et al.l ll2Q06l): (7) 
iGoicoechea&^emicha^ (12002); (8) Melnick et alj ( l201(]h ; (9) Goicoechea et alj ( l2006h : f 10) IWright et al.l boool) : (11) lEmprechtinger et al.l j2013h ; (12) IVisser et ahl ( l2013h : 113) iHerpin et al’l 
(l2012h : tl4) lBusauet et aLl(l2014ll 
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4.2.1. H 2 O abundances 


We find high abundances of H 2 O in Ccore - lower in Couter- 
In the compact inner component the H 2 O abundance is simi¬ 
lar to those found in the core component of NGC 4418 an d 
the western nucleus of Arp 220 (iGonzalez-Alfonso et aLll^oT^ . 
There are no regions in the Galaxy where an H 2 O abundance 
this high has been found over comparable scales. On smaller 
scales abundances comparable to, and higher than, that in 
Ccnre. have been found in hot cores and shocked r egions (e.g.. 


Melnick et al. 120101: lEmnrechtinger et al.l 120131: IVisser et"^ 
2013f Busauet et al. l2014t) . Comparing it to more extreme ob¬ 


jects further out in the Universe we see that Zw 049.057 has a 
H 2 O abundance exc eeding that of the H 2 O rich z = 3.91 quasar 
APM 08279-H5255 dBradford et al.ll201 ih . 

In the outer component, Couter, we find H 2 O abundances 
on the order of 10~^. This is sim ilar t o the abundances found 
by Gonzalez-Alfonso et alJ (l2010t) and iGonzalez-Alfonso et al.l 
(I2OI2I) in the warm component of Mrk 231 and the eastern nu¬ 
cleus of Arp 220, respectively. In the Galaxy, similar H 2 O abun¬ 
dances ha ve been found, for exam ple, in the denser regions of 
Sgr B2 bv ICernic haro et aP (l2006ll an d the warm shocked region 
of LI 157-Bl bv iBusauet et al.l (12'^ . 


We also find that the abundance of OH seems to vary less 
than the H 2 O abundance between Ccore and Couter- For Galac¬ 
tic sources we found fewer estimates of OH abundances, but 
there are regions around Sgr B2 and the Orion nebula that have 
abundances that are comparable to those in Ccore and Couter 
dGoicoechea & Cernicharoll2002HGoicoechea et alJl200^ . 

4.3. A Compton-thick nucleus in Zw 049.057? 

A source is, by definition, Compton-thick if the column density 
of the obscuring matter is equal to or greater than the inverse of 
the Thomson cross-section, i.e., if N// > cr:^^ ^ 1.5 x 10^"^ cm“^. 
The column density that we find in Ccore is thus enough to make 
the nucleus Compton-thick and a buried X-ray source would be 
difficult to detect. In order to penetrate the obscuring layers and 
probe the nuclear energy source in Zw 049.057 observations in 
hard X-rays or even soft gamma rays are therefore needed. 

The highly obscured nature of Zw 049.057 makes it difficult 
to unambiguously determine the source of its nuclear luminosity. 
Below we investigate two possible scenarios. 


4.3.1. Obscured AGN 


4.2.2. Why is there so much water in Zw 049.057? 


The high H 2 O abundance in the core component is a ro¬ 
bust result so we will focus our discussion here. There are 
three main scenarios which can explain H 2 O abundances this 
high: ion-neutral rea ctions starting with cos mic-ray (or X- 
ray) ionization of H 2 dyan Dishoeck et al.ll2013h. neutral-neutral 
reactions in w ^m (T > 250 K) gas 1 Nereid et al.l Il995t 
ICeccarelli et al.l[T99(x [ yan Dishoeck et alll2013h . and sublima¬ 


tion of H 2 O from dust grains as the dust temperature rises above 
~ 100 K dSandford & Allamandolal [19^ iFraser et al.1 12 OOII: 
Ivan Dishoeck et ah 2013h . If the high H 2 O abundance is due to 
the first of these processes it can persist as long as the temper¬ 
ature of the gas remains high enough. If it is due to sublima¬ 
tion from dust grains however, the abundance will eventually 
decrease unless there is an efficient formation process to counter 
H 2 O destruction by, for example, ion-neutral reactions. 

The model of Ccore indicates that the dust temperature is 
slightly higher than that required for rapid sublimation to oc- 
cur while the gas tern perature (189 + 57 K) as estimated by 
iMangum et al.l d^l3bl) is somewhat lower than that needed for 
neutral-neutral reactions to be effective. Also, we know that an 
additional, hotter, component is needed to model the continuum 
at wavelengths <10 pm. Formation of H 2 O via ion-neutral re¬ 
actions may contribute to the high abundances, but the branch- 
ing ratios in the dissociative recombination of H^O^ fav or OH 
(iHerd et alJll990l : ljensen et alj|2000l : iNeufeld et al.ll20()^ while 
we infer more H 2 O than OH in both components. 

In view of the uncertainties above, we are unable to deter¬ 
mine which one of the processes is responsible for the high H 2 O 
abundance, and a combination is most likely. One way to test the 
grain chemistry would be to search for CH 3 OH, which seems to 
be m ainly formed through hydrogenation of CO on grain sur¬ 
faces (IWirstrom et alJlTQllh . 


4.2.3. OH abundances 

The OH abundance found in Ccore is, like the H 2 O abun¬ 
dance, similar to those found i n the core of NGC 4418 and th e 
western nucleus of Arp 220 bv iGonzalez-Alfonso et al.l (l2012h . 


A possible AGN in Zw 049.057 was indicated by 
iBaan & Klocknerl (l2006l) based on radio observations and 
their nuclear classification code. With the extreme column 
density indicated by our model a hidden AGN is also allowed 
despite the relative ly weak X-ray luminosity observed by 
iLehmer et al.l (l2010h with Chandra. A way to keep a high 
co vering factor for th e obscuring matter has been proposed 
by iFabian et al.l (ll998h . They suggested a model in which a 
nuclear starburst both feeds the central black hole and hides it 
by injecting the energy needed to keep the obscuring matter in 
a highly turbulent and space-covering state. In this scenario the 
nucleus stays in the obscured state until the accreting black hole 
has grown to a certain limit were the accretion rate approaches 
the Eddington limit and the obscuring matter is blown away 
from the nucleus. 

If all of the luminosity in Ccore comes from an accreting black 
hole we can derive a lower limit to its mass of ~ 3.8 - 5.6 x 10® 
Mq by assuming that it radiates at its Eddington limit. In reality, 
due to the high optical depth of the dust continuu m the limit on 
the m ass can actually be up to 500 times higher (lEabian et al.l 
1200 ^ . 

The similarity with the western nucleus of Arp 220, sug¬ 
gested by some authors to harbor a hidden AGN, warrants 
a detailed comparison. A central concentration of hard X-ray 
emission in th e west ern nucleus of Arp 220 was found by 
[elements et al.l d2002h using Chandra X-ray imaging. Several 
possible scenarios for this emission are offered, including a weak 
AGN contributing less than 1% to the bolomettic luminosity of 
the galaxy. They also state that they cannot rule out the possi¬ 
bility that they are just seeing a small fraction of the emission 
from an obscur ed AGN and refer to the hard X-ray BeppoSAX 
observations bv llwasawa et al.l d 200 lh who found that a column 
higher than ~ 10 ^® cm“^ is required if a large fraction of the 
bolometric luminosity comes from an AGN. Based partly on the 
very high surface brightness of 5 x 10*"^ Lq kpc“^ on a scale 
of 35 pc found in their millimeter interferometer observation, 
iDownes & Eckard d20Q7h concluded that the source of the lu¬ 
minosity in the western nucleus of Arp 220 can only be a black 
hole accretion disk. They also estimated that the total proton col¬ 
umn density towards the western nucleus is ~ 1.3 x 10^® cm“^. 
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enough to hide the hard X-rays f r om an energetically important 
AGN. In addition, iTakano et al.l (l2005l) found absorption lines 
of NH 3 with full widths at zero intensit y of up to 1800 km s~ * 
possibly due to rotating gas in an AGN. ISakamoto et al.l (l2008h . 
however, do not exclude that the luminosity is primarily pow¬ 
ered by a young starburst equivalent to hundreds of super star 
clusters. 

Zw 049.057, on the other hand, w as undetected i n har d 
X-rays in the Chandra observations by iLehmer et alJ (l2010ll . 
However, as already stated, our models indicate that the col¬ 
umn density towards the center of Zw 049.057 may be as high 
as 10^^ cm~^, c ompara ble to that in Arp 220. Observations by 
iMangum et al.l (120131)1) reveals that Zw 049.057 also has NH 3 
absorption lines. These are however not nearly as wide as those 
seen in Arp 220. High spatial resolution observat i ons at millime¬ 
ter wavelengths, like those of iDownes & Eck^ (l2007h . has not 
been carried out for Zw 049.057. Doing this would be one way 
to truly test how much it actually resembles the western nucleus 
of Arp 220, and possibly also to determine if it is powered by an 
AGN or a central starburst. 


4.3.2. Starburst 


Another possible candidate for the nuc lear power source is a 
compact starburst. IVeilleux et akl (Il995h classified Zw 049.057 
as an Hit galaxy based on optical emission line ratios in 
the nucleus of the galaxy. Based o n mid-IR observations 
(IStierwalt et al.ll2013Ulnami et al.ll2013h it is also located in the 
starburst region in the diagnostic diagram bv ISnoon et al.l(l2007h 
as we ll as in mid-IR emission line diagrams (iMelendez et alJ 
12014 . 


The core component has a luminosity of 0.7 - 1.2 x 10" Lq 
and a radius of 12 - 32 pc. To put these numbers into con¬ 
text we_can_u^e_^rB2nVI) with its luminosity o f 6.3 x 10^ 
L( 7 i (Goldsmith et ^Il992h and radius of ~ 0.6 pc dRolffs et al.l 
EolF as a template. To generate the luminosity of Ccore, ~ 
(1 - 2) X 10"^ Sgr B2(M)-like sources would be required to be 
packed into a volume between 0.13 and 1.25 times that of Ccore- 


The surface brightness of the core component is very high, 
0.4 - 1.5 X 10 '"^ Lq kpc“^. This is up to one order of magni¬ 
tude higher than the limitin g surface brightness in nuclear star- 
bursts around AGN seen by IDavies et alJ (l2007h . For dust tem¬ 
peratures below 200 K this apparent limit for starbursts can be 
explained by a theoretical model in which stellar radiation pres- 
sure provides the maio rity of pressure support in the gas disk 
dThompson et al.|[2005l) . For hot starburst (Td > 200 K) this limit 
does not apply and the surface brightness ca n be higher than ~ 
10 '^ Lq kpc“^ dAndrews & Thompsonll20l'lh . The core compo¬ 
nent is not hot enough for this to happen, but our models do not 
reproduce the continuum at wavelengths <10 pm and it is pos¬ 
sible that a hot nuclear component, whic h we do not model, has 
a dust temperature > 200 K. In any case. I Andrews & Thompson! 
d 201 Ih further state that the conditions necessary to enter this 
regime may only be attained in the parsec-scale star formation 
surrounding AGN. 


Taking only the gas mass into account, our model for Ccore 
implies a luminosity-to-mass ratio of ~ 200 - 2000 Lq 
largely compatible with the limit of 500-1000 L© Mg* f or a dust- 
embedded starburst support ed by radiation pressure dScovillel 
l2003t Pfhompson et alj|2005l) . 


4.4. In- and outflowing gas 


We find a strong inverse P-Cygni profile in the [O I] 63 pm line, 
indicating inflowing gas. No obvious evidence for such inflow is 
however seen in the molecular lines so the infall might take place 
on scales larger than the ~ 100 pc molecular structure. There is 
however a possibility that the H 2 O 221 —> lio line is actually 
tracing both infalling and outflowing gas, see the discussion in 
Sect. 12.1.11 The estimated 60 micron source size is ~ 100 pc, 
so the red-shifted gas motions have to occur on these (projected) 
scales. 

A possible source of infall could be tidal debris or inter¬ 
actions where a galaxy accretes gas from a companion (as is 
the case for NGC4418, for example; Klockner, priv. comm.). 
There are several galaxies in the same held as Zw 049.057, but 
they seem to belong to a background cluster and are not com¬ 
panions of Zw 049.057. Another possibility is gas streaming in 
from a cooling flow simi lar to those seen in elliptical galaxies 
(iBregman & Athe^l2004 . Bar instabilities that acts to drive gas 
into the nucleus is a possible scenario but it seems unlikely since 
it would require rather a precise geometry. 

We also And a normal P-Cygni profile in the H 2 O 3 i 3 ^ 2 o 2 
line, indicative of outflowing gas. We do not see any outflow 
signatures in the other molecular lines however, and if the fea¬ 
ture is real it probably traces gas outflowing on the far side 
of the galaxy as discussed in Sect. 12.1.11 Compared to the 
fastest molecular outflows seen in Herschel studies of LIRGs 


(e.g..lFischer et al.ll201()t ^rm et al.ll201 lUVeilleux et al.ll2013l: 


iGonzalez-Alfonso et al.ll^l4bh the potential outflow velocities 

here are moderate, with the redshifted velocity component ex¬ 
tending to ~ 500 km s“*. 


4.5. The *'’0/'*0 isotopic ratio 

is a so called secondary nuclide, produced by He burn¬ 
ing (lMaederlll983h in massive stars and may enrich the inter¬ 
stellar medium (ISM) by winds and/or supernova ejecta. Hence 
the isotopic abundance ratio may reflect the stellar 

initial mass function (IMF) as well as the age of a burst of 
star formation. Low abundance ratios are taken to in¬ 

dicate the enrichment of the ISM by massive stars in a star- 
burst. iGonzalez-Alfonso et al.l (12012h use this ratio to place 
LIRGs/ULIRGs in an evolutionary scenario where the dust en¬ 
shrouded LIRG NGC 4418 (with an estimated ratio 

of 500) is placed in an early evolutionary phase, the ULIRG 
Arp 220 (with an estimated ratio of 70) is more 

evolved according to this scheme. The ULI RG quasar Mrk 231 
has an extreme ratio of < 30 (iFischer et al.l 1201 Ol : 

iGonzalez-Alfonso et alJl^l4bh and hence has the most evolved 
starburst of the three. A caveat here is that it is assumed that all 
three galaxies have starbursts of a similar IMF with the same im¬ 
pact and ISM mixing. Note that a large ratio may also be 

a signature of a very evolved starburst where the has caught 
up with the enrichment. 

We estimate a abundance in Zw 049.057 that is sim¬ 

ilar to or slightly higher than that in Arp 220, and this would 
therefore suggest that its starburst is in an evolutionary stage be¬ 
tween NGC 4418 and Arp 220. This is difficult to validate since 
the nuclei are so deeply obscured that standard tracers of star- 
burst evolution do not work. A possibility is to use the q-factor 
(i.e., the ratio between radio and far-IR intensity dHelou et al.l 
Il985h l - but this is aggravated by the fact that Mrk 231 is a 
quasar so its radio emission is dominated by the jet. However, 
if we take the ratios as tracing starburst evolution and 
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compare the ’’younger” NGC 4418 with the ’’older” Zw 049.057 
we find interesting similarities between them. Both have strong 
infall signatures in the [O I] 63 /im line suggesting that the cen¬ 
tral region of the galaxies are being replenished with new mate¬ 
rial that can feed the activity. NGC 4418 is interacting while this 
is not clear for Zw 049.057. Both galaxies are early type spirals 
where Zw 049.057 is showing evidence also of star formation on 
somewhat larger scales than NGC 4418. No outflow signatures 
are found in the g as of NGC 4418 (even though there are some 
hints in the dust (ISakamoto et al.ll2013h l while Zw 049.057 is 
showing tentative outflow signatures. 


5. Conclusions 

The following results were obtained from the observations of 
Zw 049.057 presented in this paper: 

- PACS spectroscopy reveals absorption in highly excited 
states of H 2 O and OH as well as in their isotopologues. 
In total, ten H 2 O lines and five OH doublets with lower level 
energies up to Slower 600 K and Slower 400 K, respectively, 
were detected in absorption. Two '^OH doublets were also 
detected in absorption. 

- The [O I] 63/rm line exhibits an inverse P-Cygni profile in¬ 
dicative of infalling gas with a velocity of 30 km s“'. A 
P-Cygni profile in the H 2 O 3i3 ^2o2 line with the redshifted 
emission extending to 500 km s ' was also detected. 

- SPIRE spectroscopy reveals a submillimeter spectrum that is 
rich in H 2 O. Eight H 2 O lines in emission and one in absorp¬ 
tion were detected along with two Hj^O lines in emission. 
The H 2 O lines had upper state energies up to Eiower 600 K. 

We have used multicomponent radiative transfer modeling 
to analyze the absorption-dominated PACS spectra together with 
the emission-dominated SPIRE spectrum and the continuum lev¬ 
els. The important conclusions from this work are: 

- Very high H 2 column density ((0.3 - 1.0) x 10^^ cm“^) to¬ 
wards the core of Zw 049.057 indicates that it is Compton- 
thick. This thick nuclear medium is also responsible for the 
frequency degradation of the intrinsic luminosity as high en¬ 
ergy photons are absorbed and re-emitted at longer wave¬ 
lengths. 

- High H 2 O column density per unit of continuum optical 
depth at 50 fim of ~ 10'* cm“^, an OH/H 2 O ratio of 0.4-0.8, 
and a dust temperature of 100-130 K are derived towards the 
same nuclear region. The molecular abundance relative to H 2 
in this region is also estimated to be high,;\fH 20 ~ 5 x 10“^. 

- The high surface brightness of (0.4 - 1.5) x lO'"' L© kpc“^ on 
a spatial scale of 10 - 30 pc is indicative of either a buried 
AGN or a very dense nuclear starburst. 

- The absorption in the '*0 isotopologue of OH indicates en¬ 
hancement of '*0 with a low '®0/'*0 ratio of 50 - 100. If 
this ratio is a tracer of starburst evolution it would place the 
galaxy in the same phase as Arp 220, in between the less 
evolved NGC 4418 and the more evolved Mrk 231. 
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